Dissipation layers in Rayleigh— Benard convection: A unifying view 
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Boundary layers play an important role in controlling convective iieat transfer. Their nature 
varies considerably between different application areas characterized by different boundary con- 
ditions, which hampers a uniform treatment. Here, we argue that, independent from boundary 
conditions, systematic dissipation measurements in Rayleigh-Benard convection capture the rele- 
vant near-wall structures. By means of direct numerical simulations with varying Prandtl numbers, 
we demonstrate that such dissipation layers share central characteristics with classical boundary lay- 
ers, but, in contrast to the latter, can be extended naturally to arbitrary boundary conditions. We 
validate our approach by explaining differences in scaling behavior observed for no-slip and stress- 
free boundaries, thus paving the way to an extension of scaling theories developed for laboratory 
convection to a broad class of natural systems. 



Buoyancy driven fluid flows arc ubiquitous in nature. 
They stir the turbulence on the solar surface power 
the zonal windpattern observed in the atmospheres of 
giant planets [H, are responsible for Earth's plate 
tectonics 0, 01 j cause atmospheric cloud formation 
and thunderstorms 0] and occur in biological systems 
[^,0. A simplified analog of such situations, the so-called 
Rayleigh-Benard configuration consisting of a plane fiuid 
layer heated from below, served as a cornerstone for the 
development of hydrodynamic stability theory [lol - [l^ 
and has become a paradigm for studies of convective tur- 
bulence. 

The key quantities for characterizing convective sys- 
tems are the resulting heat and momentum transport, 
typically expressed in terms of the non-dimensional Nus- 
selt number Nu and Reynolds number Re. Much re- 
search over the past years has focussed on the prediction 
of these transport characteristics, with the Grossmann- 
Lohse (GL) scaling theory [l^ and its extensions [l^ 
[lij being among the most prominent examples. Beyond 
their significance in testing our theoretical understand- 
ing against experimental evidence, such scaling theories 
are crucial in estimating the heat transport in natural 
situations, for which the control parameter values often 
differ from those accessible in laboratory or numerical 
experiments by many orders of magnitude [l7j . 

Modern scaling theories 13 , emphasize the impor- 
tance of both the thermal and viscous boundary layer 
in controlling the scaling behavior. While revealing 
an irnpressive consistency with available experimental 
data [l9| , these theories unfortunately cannot be applied 
directly to natural situations like the ones mentioned 
above. The problem is that the classical picture of a 
viscous boundary layer, in which the tangential veloc- 
ity components decrease rapidly over a smaU, 0(Re^/^) 
length scale toward the boundary, is inextricably linked 
to the presence of a rigid boundary surface. In most of 
the examples mentioned above, however, there is no rigid 
boundary forcing the tangential velocities to drop to zero. 



and consequently no classical viscous boundary layer is 
to be expected. Theories relying on the presence of such 
a layer, including the scaling theories mentioned earlier, 
therefore cannot be generalized to arbitrary boundary 
conditions in a straightforward manner. 

In this Letter we introduce the concept of what we call 
dissipation layers, a generalization of the classical bound- 
ary layers that is based on the kinetic energy budget 
and the thermal variance balance. To test this concept, 
we study these dissipation layers for both, laboratory- 
style no-slip boundary conditions and so-called stress- 
free conditions, in which the horizontal shear stresses are 
required to vanish at the boundary. The classical bound- 
ary layer picture only applies to the former, but not to the 
latter case, in which the horizontal velocities tend to peak 
on the boundary itself. Interestingly, we find that pro- 
nounced dissipation layers are observed for both types of 
boundary conditions and in many respects behave qual- 
itatively similar. This paves the way for applications of 
the existing theories to natural systems. 

Since the boundary layer dynamics is controlled by 
molecular diffusion processes, it is very sensitive to 
changes in Prandtl number Pr = i^/k, where k and v 
are the thermal and momentum diffusivities. Pr varies 
strongly for different fluids, ranging from Pr — O(10~^) 
in stellar plasmas to Pr = 0(10^"^) in planetary mantles. 
Both of the above reasons have motivated several nu- 
merical studies (iol-lis} on the influence of Pr on convec- 
tion. Surprisingly, for stress-free boundary conditions, 
scaling laws similar to the ones proposed by Grossmann 
and Lohse have been observed [24I . Pointing out the ab- 
sence of classical viscous boundary layers for stress-free 
boundaries, previous authors [23] thus concluded that 
these do not control the observed scaling and thus casted 
doubt on the universal applicability of the GL-theory. 

A number of important questions immediately arise in 
this context: i) Is there, for stress-free conditions, a vis- 
cous layer close to the boundary [H, US] that plays a 
role similar to the classical viscous boundary layer in the 
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FIG. 1. Temporally averaged depth-profiles of the horizon- 
tally averaged viscous dissipation rate (red solid line) and the 
thermal dissipation rate (blue dashed line) for no-slip (upper 
graph) and stress-free (lower graph) boundaries and for dif- 
ferent Prandtl numbers. The horizontal axis is scaled by the 
corresponding globally averaged dissipation rates (cf. equa- 
tion Q). The thickness of the newly defined dissipation layers 
are illustrated for each Prandtl number and boundary condi- 
tion, respectively. 



nO"Slip case? ii) If so, is there a unifying way of defin- 
ing near-wall viscous layers independently of the exact 
nature of the mechanical boundary conditions? iii) Can 
this unified definition also be extended to describe the 
thermal boundary layers? iv) If this is the case, are the 
properties of these newly defined layers compatible with 
the basic assumptions of the established theories? v) Fi- 
nally, can these layers be connected to specific scaling 
regimes for Nu and Re, effectively allowing a generaliza- 
tion of the existing theories? 

Here, we take the view that the most basic feature of a 
boundary layer region is that molecular dissipation pro- 
cesses play a major role there, as has been demonstrated 
in previous work for laboratory style, no-slip boundaries 
[25I [2^ . We thus focus on the horizontally averaged ki- 
netic energy dissipation rate. 



(6„)Jz) = -2Pr(Tr(S2)), , 



(1) 



defined by the trace of the squared rate-of-strain tensor 
S, and the corresponding thermal dissipation rate. 



(eT),W = (-(VTf). 



(2) 



Using these quantities, the edges of the viscous and ther- 
mal dissipation layers arc defined by the vertical posi- 
tions where the local viscous and thermal dissipation 
rates equal their volume-averaged values, 

{^u)h (^u,dl) = and {eT)^ (^t,dl) (eT)y • (3) 
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FIG. 2. Thickness of the dissipation layers versus the Prandtl 
number for no-slip (upper graph) and stress-free (lower 
graph) boundaries. Red circles represent the thermal dissi- 
pation layer while blue squares denote the viscous dissipation 
layer. The gray arrows indicate the crossover of the thick- 
nesses of both dissipation layers (At,dl = Au,dl). The classi- 
cal thickness definitions of the thermal and viscous boundary 
layer, At,bl and A„,bl, are denoted by green triangles and 
orange crosses, respectively. Note that for stress-free bound- 
aries, the classical viscous boundary layer definition cannot 
be applied. 



The distances from the closest boundary then define the 
viscous and thermal dissipation layer thicknesses A„_dl 
and At,dl. 

Note that this definition, which separates the flow into 
regions of high (above average) and low (below average) 
thermal and kinetic energy dissipation, is independent 
of the nature of boundary conditions. We also focus on 
dissipation rates here because these are expected to be 
related to the overall heat transport; in particular we 
have [ill 



(e«)y = (Nu - l)RaPr and (eT)y = Nu 



(4) 



In the following, we investigate the properties of these 
dissipation layers in detail. Three-dimensional direct nu- 
merical simulations of Rayleigh-Bcnard convection in a 
Boussinesq fluid with periodic boundary conditions in the 
horizontal direction have been carried out using an accu- 
rate pseudo-spectral method [27[. The top and bottom 
boundaries are impermeable, kept at a fixed tempera- 
ture, and either the horizontal velocities (no-slip) or the 
shear stress (stress-free) are assumed to vanish at the 
boundary. For a fixed Rayleigh number Ra = 5 x 10^, 
the Prandtl number range 0.01 < Pr < 300 is systemati- 
cally explored. Even for this moderate Rayleigh number, 
spatial resolutions up to 576'^ grid points were necessary 
to adequately resolve the Kolmogorov scales within the 
bulk at low Pr [ii,!!^. 

Figure [1] shows the height-resolved dissipation rates for 
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no-slip and stress-free boundaries and varying Prandtl 
numbers. All cases exhibit regions of strongly enhanced 
dissipation close to the boundary, validating the concept 
of dissipation layers separating the boundaries from the 
bulk. The profiles and consequently the layer thicknesses 
strongly depend on Pr: While for low Prandtl num- 
bers the viscous layer is much smaller than the thermal 
layer, the opposite is observed for high Prandtl num- 
bers. A crossover of the layer thicknesses takes place 
around Pr = 1 for no-slip and slightly below for stress- 
free boundaries (cf. also FIG. [2]). This "change of hierar- 
chies" separating the parameter space into two distinct 
regimes has been predicted for the no-slip case in previ- 
ous theoretical works [3, ■ 

The dissipation layer thicknesses along with their clas- 
sical counterparts At,bl and A„.bl are shown as a func- 
tion of Pr in FIG. [5] The classical boundary layer thick- 
nesses are defined using linear fits of the depth profiles 
in the vicinity of the boundary, i.e. the so-called slope 
method [s^l- It is observed that, apart from a con- 
stant prefactor in amplitude, the viscous dissipation layer 
shows a similar Prandtl number dependence as the classi- 
cal viscous boundary layer in the no-slip case: The layer 
thickness increases with Pr and starts to saturate for the 
highest Prandtl numbers. The thermal dissipation layer 
and classical thermal boundary layer also show similar 
behavior and decline slowly with Pr. Differences appear 
for the lowest Prandtl number, which can be explained 
by the fact that the thermal boundary layer in this case 
is less well defined due to a finite thermal mean gradi- 
ent in the bulk. We conclude that the newly defined 
viscous and thermal dissipation layers capture the near- 
wall characteristics of the system equally well as the clas- 
sical approaches, which however fail to provide a concise 
definition for a viscous boundary layer in the stress-free 
case. 

We now demonstrate that the new framework indeed 
allows a generalization to stress-free boundary condi- 
tions. As a central result we observe in FIG. [2] that the 
layer thicknesses show a functional form very similar to 
the no-slip case. The viscous dissipation layer thickness 
increases with Pr with a tendency to saturate at high 
values, whereas the thermal dissipation layer thickness 
decreases with Prandtl number consistently with its clas- 
sical counterpart. Again, the change of dissipation layer 
hierarchies is clearly visible for stress-free boundary con- 
ditions. 

The usual way to estimate the viscous boundary layer 
thickness is to apply the classical Prandtl-Blasius the- 
ory to convective systems resulting in the pre- 

-1/2 



diction A„ ~ Re 



34|. 



To check whether or not this 
prediction also holds for the proposed dissipation based 
thickness definitions, FIG. [3] shows Au^dl as a function 
of Re for both stress-free and no-slip boundary con- 
ditions. For convenience, the classical viscous bound- 
ary layer thickness A„,bl is also shown for the no-slip 
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FIG. 3. Thickness of different definitions of a viscous layer 
versus tlie Reynolds number Re = a/ {u'^)v h/u, where h is 
the height of the fiuid layer and u the momentum diffusivity. 
The orange crosses denote a classical definition of a viscous 
boundary layer thickness Au,bl, while the blue squares and 
the red circles denote the dissipation layer thickness A„,dl for 
stress-free and no-slip boundary conditions, respectively. 
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FIG. 4. Ratio of the DL contribution of the globally aver- 
aged kinetic dissipation rate tu versus the Prandtl number. 
The blue squares denote results of simulations with stress- 
free boundaries, while red circles denote results with no-slip 
boundary conditions. The gray arrows divide the Prandtl 
number range into a bulk and a dissipation layer dominated 
regime. For completeness, the decomposition employing the 
classical boundary layer definition is represented by orange 
crosses. 



case. Fits to the numerical data for large Re yields 

A„,BL ~ Re-°-4'^6±0.004 ^^^^^ ^ j^g-0.474±0.011 

the no-slip case, both in fair agreement with the classical 
Prandtl-Blasius prediction. In contrast, for the stress- 
free case, A„ dl decreases faster with Re, resulting in a 
scaling law A„,DL ~ Re-° *^"" "^^ 

Within the parameter range covered by our study the 
system shows a transition from bulk dominated viscous 
dissipation at low Prandtl number to a regime charac- 
terized by significant dissipation within the dissipation 
layer at high Prandtl number. Figure |4] shows the ra- 
tio of the dissipation layer contribution (e„)DL to the 
globally averaged kinetic dissipation rate {eu)v versus 
the Prandtl number. In the low Prandtl number regime. 
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(eu)DL/(eu)y is small and therefore a strong dominance 
of the bulk contribution is observed for stress-free and 
no-slip boundaries, respectively. The dissipation layer 
contributions then increase with Prandtl number, before 
saturating at Pr = 0(10) with a contribution of approx- 
imately 50% for the no-slip case and 35% for the stress- 
free case. The saturation point marks a clear change in 
dynamical behavior, and, borrowing from the terminol- 
ogy of Grossman and Lohse, we use it to quantify the 
transition from bulk to dissipation layer dominance. 

On the basis of our newly defined dissipation layers, 
we can identify three different regimes within the inves- 
tigated Prandtl number range, with the limits being de- 
fined by the crossover of the thermal and the viscous 
dissipation layer thicknesses and by the saturation of the 
dissipation layer contribution to {tu)v other words, 
the regimes are classified by the "hierarchy" of the vis- 
cous and thermal dissipation layers and by the relative 
strength of viscous bulk and near-wall dissipation. To 
illustrate that these regime transitions have dynamical 
relevance, FIG. [5] shows the Nusselt number obtained 
from numerical experiments versus the Prandtl number. 
In the stress-free case, we find three clearly distinguish- 
able scaling regimes, each characterized by a different 
Nusselt-Prandtl scaling law. Remarkably, the observed 
transitions between these scaling laws agree well with 
the regime limits described above. For no-slip bound- 
aries, the intermediate Prandtl number regime appears 
to be very narrow. It is unclear whether the correspond- 
ing data points reflect a distinct power law, or merely 
represent a gradual transition between high- and low 
Prandtl number scalings. This observation is consistent 
with the fact that our regime classification also suggests 
a very narrow intermediate Prandtl number regime. It is 
also consistent with the theory of Grossmann and Lohse, 
which, for the Rayleigh numbers used in this study, pre- 
dicts an intermediate regime so narrow that the authors 
are led to question its very existence in their original 
work [l3| . Future simulations at higher values of Ra are 



Cji bulk dominates C 



needed in order to obtain more conclusive data. 

The results presented in this Letter answer many of the 
questions posed at the beginning of this paper. In partic- 
ular, we have shown that near walls, narrow regions of en- 
hanced viscous and thermal dissipation are a generic fea- 
ture of turbulent convection irrespective of the exact na- 
ture of the mechanical boundary condition. This observa- 
tion allows for a definition of thermal and viscous "dis- 
sipation layers" that neither relies on the applicability, 
nor on the validity of classical boundary layer models. In 
cases where these dissipation layers have classical bound- 
ary layer counterparts, both layers have been shown to 
exhibit a similar Prandtl number scaling, albeit with a 
different prefactor. Different from their classical coun- 
terparts, however, the thermal and viscous dissipation 
layer thicknesses are shown to cross around Pr = 0(1) 
for both, no-slip and stress-free boundary conditions. 
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FIG. 5. Nusselt number versus the Prandtl number for no- 
slip (upper graph) and stress-free (lower graph) boundaries. 
The vertical dashed line indicates the "change of hierarchies" , 
i.e., the Prandtl number where the thermal and the viscous 
dissipation layers are of the same thickness. The color coded 
background represents the ratio of the DL contribution of e„ 
in FIG. 131 where light gray indicates a bulk dominated regime 
and dark gray indicates a dissipation layer dominated regime. 
Data fits for the no-slip case are in good agreement with the 
predictions by Grossmann and Lohse [31 (^ Pr^^^ for low 
Prandtl numbers) 



This separates the parameter space into distinct regions 
characterized by different dissipation layer hierarchies, in 
accordance with the theoretical predictions [l3, 3 29j . 
The contribution of the dissipation layers to the overall 
dissipation has also been shown to be consistent with the 
theoretical assumptions. Finally, a regime classification 
based on dissipation layer hierarchies and on their con- 
tribution to the overall dissipation correlates well with 
transitions observed in the Nu — Pr scaling. This strongly 
suggests that existing scaling theories developed for lab- 
oratory convection can indeed be extended to boundary 
conditions relevant for a broad class of natural systems. 

An interesting finding in this context is that for stress- 
free boundary conditions, the viscous dissipation layer 
thickness decreases more rapidly with Re than expected 
from a balance between horizontal advection and verti- 
cal diffusion, which would lead to the classical A„^dl ^ 

I/O 

Re ' scaling. Although an explanation of this obser- 
vation currently remains elusive, it might have strong 
implications for scaling transitions caused by the change 
in dissipation layer hierarchy. We might speculate that 
future studies of convection with stress-free boundaries 
are likely to reveal exciting surprises, such as transition 
scalings different from those predicted for no-slip systems 
by the current theories 35|, |36[ . 
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